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Entropy Production in Gluodynamics 
in temporal axial gauge in 2+1 dimensions 

Akihiro Nishiyama and Akira Ohnishi 
Yukawa Institute for Theoretical Physics, Kyoto University, Kyoto 606-8502, Japan 

Entropy production in non-Abelian gauge theory is discussed in the limit of vanishing 
classical fields. Based on the Kadanoff-Baym equation with the leading order self-energy in 
the temporal axial gauge, the kinetic entropy is introduced and is proven to increase for a 
given self-energy. Numerical simulation of the Kadanoff-Baym equation for the non-Abelian 
gauge theory is performed in 2+f dimensions. Starting from anisotropic distribution in 
momentum space, relaxation to isotropic distribution proceeds and the entropy is produced 
due to off-shell scattering of massive transverse polarization mode as expected from the proof 
of the H-theorem. 

& ' 81. Introduction 

PLh ■ 

In heavy-ion collisions at the Relativistic Heavy Ion Collider (RHIC), it is plau- 
sible from theoretical arguments and experimental evidences that a new form of 
matter made of quarks and gluons is created.® This matter referred to as a strongly 
coupled QGP (sQGP)® is characterized by its strongly interacting nature; the shear 
viscosity is extremely small and the thermalization time is very short. Hydrody- 
namic models successfully describe the radial and elliptic flows at RHIC®'® and 
they require a local thermal equilibrium in the early stage. Comparison with the 
experimental results of the elliptic flow suggests that a thermalization time of the 
■ partons is about 0.6-1 fm/c.® This short thermalization time seems to contradict 

the perturbative estimate,® and it is compatible with the formation time of partons. 
This means that the Boltzmann equation would have serious problems in describing 
t-» the thermalization of dense systems, such as glasma, formed in the early stage of 

collisions. 

At lower energies, hadronic transport models based on the Boltzmann equation 
have been successfully applied to heavy-ion collisions including nuclear fragment 
formation,® flavor production,® particle spectra,® and collective flows.® By com- 
parison, hadronic transport models fail to describe early thermalization and predict 
smaller elliptic flows at RHICP® These successes and failures imply the formation 
of QGP at RHIC, and may also suggest the limitation of the Boltzmann equation in 
dense systems. Boltzmann approaches are based on the assumption that the parti- 
cle interaction is rare and three particle collisions can be ignored, and the spectral 
function is narrow enough for the quasi-particle approximation to be valid. In dense 
matter formed at RHIC, one or both of these conditions are not satisfied. One of 
the prescriptions to implement multi-particle interaction is to apply the classical 
field theoriesP®"^® "When the gluon density is so high that gluons are combined 
to form finite color glass condensate, a large part of the energy density is stored 
in the classical field. It is generally believed that the early thermalization is real- 
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ized from rapid growth of intense classical Yang-Mills fields triggered by the plasma 
inst abilities IID-tlDJISJilIS The decay of strong gluon field into quarks and gluons 
needs quantum field description in nonequilibrium conditions. In addition, classical 
field dynamics should not be applied to high momentum modes even in the dense 
early stage. Another aspect of dense systems, the spectral function change, can be 
described perturbatively or by using some kind of resummation in quantum field the- 
ories. Among various ways of resummation, the ^-derivable approach would be the 
most favorable, since it does not have the secularity problem, and leads to a correct 
equilibrium. Hence we should study the gluon thermalization processes based on 
nonequilibrium quantum field theories such as the Kadanoff-Baym equation based 
on ^-derivable approximation. 

The Kadanoff-Baym (KB) equation, equivalently the Schwinger-Dyson equation, 
is given by Kadanoff and Baym®'® i n the reformulation of a functional technique 
with Green's functions by Luttinger and Ward.® It is related to variational tech- 
nique of an effective action using the so-called ^-derivable approximation, and is 
shown to satisfy the conservation laws of a system.® In the ^-derivable approxi- 
mation, diagrams contributing to the self-energy are truncated with some expansion 
parameter. This technique is applied to relativistic systems and reformulated for 
nonequilibrium many body systems by using the Schwinger-Keldysh path integral 
formalism and 2-particle irreducible (2PI) effective action technique.® - ® 

Nonequilibrium quantum field theoretical approach is a tool to understand a 
large variety of contemporary problems in high energy particle physics, nuclear 
physics, astrophysics, cosmology, as well as condensed matter physics.®'® In cos- 
mology it is applied to describe the reheating process due to the decay of inflaton 
field coupled to the Higgs at the end stage of inflation.®"® In condensed matter 
physics it is applied to Bose-Einstein Condensate. In this approach both conserva- 
tion law and gapless Nambu-Goldstone excitation® are respected.®'® It might 
be also possible to apply it to heavy ion collisions to describe equilibration processes 
of par tons. 

With respect to numerical simulation, Danielewicz first carried out a pioneering 
work in 198<P^ to describe heavy ion collisions with spectral functions at the non- 
relativistic energies. In the relativistic level, the equilibration processes are simulated 
in (p 4 model in the symmetric phase in 1+1®® 2+lE2} anc i 3+iHQ},IIiJ di mens i ns. 
Numerical analyses have been also done in the broken phase in 1+1® and 3+1® 
dimensions. Recently simulations have been extended to the system in an expanding 
background.® Similarly based on the next-to- leading order self-energy in the sys- 
tematic 1 /N expansion,® simulations have been done both in the symmetric® and 
broken® phases in the O(N) model. All of these analyses show that the thermal- 
ization proceeds; distribution functions lose information on the initial condition and 
converge to the Bose-Einstein distribution. Until now the KB equation for the non- 
Abelian gauge theories has never been solved numerically. This is because we have 
several problems such as gauge covariance of the equation of motion,®'® infrared 
singularity, and violation of Ward identity.® 

In this paper we discuss the gluon thermalization in the temporal axial gauge 
(TAG) on the basis of the KB equation with the 2PI effective action in the lead- 
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ing order of the coupling expansion. We expect that the off-shell particle number 
changing processes such as g f-> gg contribute to the entropy production in heavy 
ion collisions. These off-shell effects come from the gluon decay due to finite spec- 
tral width and memory integral in the KB equation. A part of gg <-> gg processes, 
99 9 ^ 99i is included in the present treatment via the sequential off-shell pro- 
cesses of gg g and g f-> gg, and contributes to thermalization. In the Boltzmann 
equation, g <H> gg is prohibited,® and the on-shell gg O ggg process is the most ef- 
ficient in entropy production. It is known that the on-shell dynamics predicts a long 
equilibration time r eq = 2 — 3fm/cP and it can not explain the early thermalization 
of glasma. Hence we investigate off-shell dynamics by using KB equation and give 
whether thermalization occurs or not due to the effects which have been neglected. 
The present entropy production mechanism from the off-shell gluons is different from 
the classical plasma instabilities, such as the WeibeP^^H}'^ and Nielsen-Olserfl^J 
types. Our approach is also different from classical statistical approach on a lattice, 
which conjectures the significance of the leading order and next-to-leading order 
resummed loop self-energyP3 The aim in this paper is to show that the off-shell 
processes g — >■ gg can contribute to equilibration of gluons. For this purpose, we 
introduce the concept of kinetic entropy and the H-theorem, and perform numerical 
simulations of the KB equation in gluodynamics. In order to avoid the problems 
mentioned above, we make some prescriptions and approximations, such as gauge 
fixing, infrared cutoff and simulation of only the transverse scattering contributions. 

This paper is organized as follows. In Sec. [21 we write down the KB equation 
in gluodynamics for both transverse and longitudinal part of Green's functions. As 
the self-energy, we take account of the leading order (LO) terms in the coupling 
expansion, which contain the off-shell g f-> gg processes. In Sec. O we introduce 
kinetic entropy, and by use of the gradient expansion as in Refs. [50|) - I53p we show 
an analytical proof of the H-theorem for the given self-energy in order to examine 
whether the above leading order off-shell processes (g -H- gg) can contribute to equi- 
libration or not. The concept of entropy should also be extended to the systems 
with strong classical fields,^ 1 which is out of the scope in this paper. In Sec. HI we 
give numerical results of the Kadanoff-Baym equations for the transverse scattering 
processes with massive mode in a spatially uniform 2 + 1 dimensional spacetime in 
order to extract the effects of the off-shell particle number changing processes. We 
find that the off-shell particle number changing processes can actually contribute to 
entropy production and thermalization. We summarize our work in Sec. EJ 

§2. Kadanoff-Baym equation for gluons in Temporal Axial Gauge 

In this section we write down the Kadanoff-Baym equation for gluons in the 
temporal axial gauge (TAG), Aq = 0. This section focuses on the practical use of 
the KB equation. Let us neglect the classical field {A) = and concentrate on the 



*' The process g «->• gg contributes in the leading order 0(g 2 ), while gg «-> gg and g <H> ggg are 
the next-to-leading order 0(g 4 ) contributions of the coupling expansion in the ^-derivable approxi- 
mation. In this paper we concentrate on the leading order effects. 
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Fig. 1. Integral path in the time coordinate. 



quantum fluctuations. We start with the Lagrangian density 



C = --F a ^F^, A a Q = 0, (2-1) 



where = d^A^ — d v A°^ + gf abc A b ^A^ and f abc is a structure constant in non- 
Abelian gauge theory. The Greek indices run over the coordinate space {0, 1, • • • , d}, 
and the color indices {a, &,•••} run over {1, • • • , N 2 — 1}. We use a closed time path 
C along the real time axis, the path from to to oo and from oo to to as shown in Fig. [H 
in order to trace nonequilibrium dynamics.^'^ Then the 2PI effective action with 
vanishing classical field ({A) = 0) is written as 

r\D] = ^Trln^r 1 + % -D^D + \r 2 [D). (2-2) 

Here iF)^\- = — (5ijd x + d x ^d x ^)5 ab 5c{x — y) is the free Green's function, and D 
(D 00 = D 0i = D i0 = 0; ■■■ = 1, ■ ■ ■ ,d) is the full Green's function, both of 
which are defined on the closed time path C. The remaining functional ^/^[-D] 
in (|2-2p is a sum of the all possible 2PI graphs written with respect to D which 
remains connected upon cutting two Green's function lines. The stationary relation 
for the effective action (|2-2p . 5r/5D = 0, gives the Schwinger-Dyson equation for 
the Green's function D{x, y) 

D- 1 (x,y) = D^ 1 (x,y)-n(x,y) (2-3) 

or equivalently, 

(S ik d 2 x + d x ^d xM )Dt b {x, y) + i f dzll%(x, z)D cb (z, y) = -i6 ab 6ijSc(x ~ v) i 2 ^) 

with the proper self-energy defined as II = i5r2[D]/5D. The self-energy is expressed 
as 

II(x, y) = -iSc(x - y)II\ oc (x) + II non i(x, y) (2-5) 



where the local part iJj oc contributes to the mass shift while the non-local part 7T non i 
induces the mode-coupling between the different wavenumbers. 
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We are now interested in {ij} components of the Green's function and the self- 
energyj^ and shall decompose the two-point function Df b (x,y) and the non-local 
part of the self-energy n® b , • ■ into the statistical and spectral parts, 
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(2-6) 
(2-7) 



The upper label "1" represents the path from to to oo and "2" represents that from 
oo to to in the closed time path contour C. The spectral function p{x,y) contains 
the information on which states are the most realized, while the statistical function 
F(x, y) determines how much a state is occupied. The Schwinger-Dyson equation 
(|2-3|) can be equivalently rewritten in terms of F°- b {x,y) and pf b (x,y) as coupled 
integro-differential equations 



(<%«9* + d x ,d xJ )5 ab + n^{x) F$(x, y) 



ab 



->6c / 



dzR 



ab 
F,ij 



j (x,z)p bc k (z,y) 



i-x" 

/ dzn% 3 (x,z)F%(z,y), (2-8) 
J t 



(Sijd* + d x>i d x>j )8 ab + n^{x) pf k (x, y) 



ab 



dzmUx,z)p%{z,y) (2- 



be , 



where io is the initial time. The set of equations (|2-8|) and (|2-9|) is the Kadanoff-Baym 
equation of the gauge theory, which describes the time evolution of the fluctuations 
with F and p. At each time step the spectral function pfj(x, y) must satisfy the 
conditions following from the commutation relations: 

- abr ~ ~v -o, 

<%<^ d (x-y), 
0. 



Pij(x,y)\x°^y° 



dx0pij(x,y)\ x o^ y 

d x od y op? b (x,y)\ x0 _, y0 =0. (2-10) 

In this paper we restrict ourselves to the spatially homogeneous and color isotropic 
situation. From the translational invariance, we make a Fourier transformation of 
the Green's functions and the self-energies with respect to the spatial relative coor- 
dinates, and decompose into the transverse and longitudinal components aj*^l 
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ab 



I T (^, y ;k) + ^/ L (x , y °;k) 



(2-11) 



where I = D,II, F, p, FIp,Fl p or i7i oc . Then the KB equations of the transverse and 
longitudinal parts are simplified in the momentum space as 



[<9 2 + k 2 + 77 loCiT (x )] F T (x°, y°; k) 



[n FtT p T ] y 



to 



[d 2 xfl + k 2 + 77 loc , T (z )] p T (x°, y°; k) = - [iT^rfc 



(2-12a) 
(2-12b) 



*•* We adopt the notation 



/2 Dl instead of l 
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[<t + 77 loc , L (x )] F L (x ,y ,k) = {[n F)L p L ]y° - [n PtL F L ]£} , (243a) 

[dl Q + n loc , L (x )] P L(x°,y°;k) = - [n p , L p L ] x y ° , (2-13b) 
where we have introduced a short-hand notation, 

[AB]ll= dz°A(x ,z ;k)B(z ,y°;k) . (2-14) 
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Fig. 2. Leading order self energy for gluons. 

We approximate the self-energy 77 with the skeleton diagrams obtained at the 
leading order in g 2 which is depicted in Figj2j The transverse and longitudinal parts 
are explicitly given in terms of Dt and Dl as 



77 



\oc,T{X 



0\ 



g 2 N 



d d p 
(2vr) d 
id 



(d - 2)D T + D L + S -^±(D T - D L ) 
d — 1 



ni oc , L (x°) = g 2 N J ^ [(d - 1)D T + (D L - D T ) sin 2 p , k ] 



(2-15) 
(2-16) 



where the arguments of D in the integral is (x°,x°; p). Similarly the transverse and 
longitudinal components of the nonlocal part 77 non i are given as 



iTnoni,r(x°,y°;-k) 



9 2 N 



2(d 

g 2 N f d d p 



IT J (xF ^ P ^(P' ~p ~ k )A*(p)£>/?(-p - k ) > 

(2-17) 



a,/3 



g ' TV f d p , 

n non i )L (x°,y°;-k) = — / — — ^^Q a/3 (p, -p - k)D a (p)Dp(-p - k) , 

(2-18) 

where the sum runs over af3 = TT, LT, TL and LL. In the integral, the temporal 
variables for D a and Dp are (x°,y°). Here P a p and Q Qj g are given as follows 



p tt (p, q) =pg( P , q) + -p - q) + Ptt(-p - q, p 



QTr(p,q) 



, [(p + q)-(p-q)] ; 
(p + q) 2 



d- 2 + 



(p-q) 

p2q2 



21 



(2-19a) 
(2-19b) 
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PLi(p,q)=(p + qr 1 
Plt (p,q) =Ptl (q,p) = 

<3LT(p,q) =QTL(q,p) = 
QLL(p,q) =o , 



where Pj^rf is given as 



P$(p,q) =4(d-l) 



P 



i (p+q)2 



i 



IP • qj 

p2q2 



Q TT (q, -p - q) 
: J°ii(P,-P - q) 



(p-(p + q)) 2n 

(p + q) 2 
(p • (p + q)) : 

p 2 (p + q) 2 



(2-19c) 

(2-19d) 
(2-19e) 
(2-19f) 



2h r 



1 



(q • (p + q)) 



21 



(2-20) 



q 2 (p + q) 2 

We find several characteristic features of these coefficients. First, P a $ and Q a p 
are all semi-positive definite for d > 2, 

P al 3 > , Q a/3 > (d > 2) . (2-21) 

Next, we find some symmetric relations. For example, Ptt is symmetric for any 
exchange of the variables p, q and k, which satisfy p + q + k = 0, 

Ptt(p, q) = P TT (q, k) = P TT (K p) = Prr(q, p) • (2-22) 
By using the identity, 

2 [p-(p + q)] 2 2 [q-(p + q)] 2 



=q 



(2-23) 



(p + q) 2 (p + q) 2 

we find that PjQ is symmetric under the exchange of p and q, and the cyclic sum 

of -Pyy leads to the cyclic symmetry. Other symmetric coefficients (Qtt an d Pll) 
are also symmetric under the exchange of two momentum, 

Qtt(p, q) = Qrr(q, p) , Pll(p, q) = P LL (<k p) • (2-24) 

For mixed coefficients (Plt, Ptl,Qlt,Qtl) Eqs. (j2- 19d|) and (|2-19e|) are rewritten 
by using these vectors (p, q, k) as, 

Plt(K p) = Ptl(p, k) = Qtt(p, q) , Qlt(p, q) = QrUq, p) = Pll(p, k) . 

(2-25) 

These relations are useful to prove the H-theorem as discussed later. 

In the limit of k — > 0, we can reproduce the mass shift derived by using the 
Schwinger-Dyson equation in thermal equilibrium,^ 

d d p 



iT ri (x°,y°) =g z N(d-l) 



(2vr) d 



i5(x» - y u )[(d- l)D T (x u ,y"; p) + D L (x^ y u ; p)] 



p 2 D T (x°, y°; p) [2D T (x°, y°; p) + D L (x°, y°; p) 



(2-26) 



where Tin = (d — + II l- This result is consistent with the hard thermal loop 

(HTL) approximation As already discussed in Ref . I47p . the Ward-Takahashi 
identity is not satisfied in the KB equation with 2PI effective action, because dia- 
grams with different orders in g 2 are resummed. 
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§3. Proof of the H-theorem for the non-Abelian gauge theory 



In this section we give the expression of the kinetic entropy in terms of the 
two-point Green's functions D(x,y) for the non-Abelian gauge theory. We shall 
investigate the case with vanishing classical fields {A) = and color isotropic Green's 
functions D ab (x,y) = S ab D(x,y) in the temporal axial gauge Aq = 0. Then we find 
that the kinetic entropy satisfies the H-theorem. To be generic, we here consider 
non-uniform system. In the leading order of the gradient expansion, we can still 
utilize the same coefficients P a p and Q a p as those shown for the uniform system 
in the previous section, since the modification of these coefficients are in the higher 
order in the gradient expansion as discussed later. 

The KB equations for the transverse and longitudinal components, Eqs. (|2-12a|) - 
(|2-13bp , have a similar form to that for the scalar field theory, then the kinetic entropy 
current would have a similar form to that derived in the scalar theoryF^^lS'^S'^HJ 
We start from a conjecture that the kinetic entropy current per color in the gauge 
theory with temporal axial gauge is given by 



d d+1 k . 



+ 



(27T) 

r d d+1 k 

J (2^ 



d+l 



1) 



_ 1 dRe n T ji 
2 dk» 
1 ORe 77 L , R 
" 2 dk u 



pr_ 1 dRe D t ,r H P ,T 
i 2 dk^ i 



Pl_ 1 dRe £>l,r I7 PiL 
i 2 dk„ i 



(3-1) 



where 



a(/) = (l + /)log|l + /|-/log|/| , 



(3-2) 



R = i(D^-D L c f),n a ^ R 
ipafa and D 2 ^ = -ip a (l + f a 



— 11^ ) and f a is defined with expressions 



12 



(a = T,L). We begin with the Schwinger-Dyson 
equation Eq. (|2-3p . In the similar way as in Ref . I52"|) , we obtain the following relations 
with respect to the Fourier transformations D(X,k) = f d(x — y) e ip ( x ~ y ^ D(x, y) 
where X 



is the " center-of-mass" coordinate. 



d d+l k 



D 



12 



1 



(27T) 



^T{d-l)C T \og^ 1+ - 



d d+1 k n 12 

^Cd^ (3-3) 



(27T) 



In deriving of the above relation we have used gradient expansion and taken the 
1st order of the approximation. (If the singularity in the longitudinal mode ~ qjjjz 



appears, it cancels in the fraction of logarithmic part.) See, for example, Ref. 
for the detailed derivation of the entropy current and its divergence, Eqs. (|3Tp and 
(|3-3p . from the KB equation in Eq. (|2-3p . Here the term C represents the collision 
terms, 



C T (X, k) 
C L (X,k) 



i[n T , P {x,k)F T {x,k) 
i[n LtP (x,k)F L (x,k) 



n T>F (x,k) PT (x,k)] 
n L>F (x,k) PL (x,k)] 



(3-4) 



The difference from the scalar theory is that the kinetic entropy is the sum of the 
transverse and the longitudinal part. 
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The remaining work is to prove that the R.H.S. of (j3-3|) is positive definite. It is 
convenient to express IIt(X, k) in the collision term of (|3-3j) . which is obtained from 
the Fourier transforming with respect to the relative coordinate x° — y° and is given 
as 



n T (X, -k) = - j J2 P ^(P> "P " *)D a (X,p)Dp(X, -p - k). 

(3-5) 



In a similar way, the longitudinal part TIl{X, k) can be expressed by changing 
D(x° ,y° ,p)D(x° ,y°; -p - k) to f -^D(X,p)D(X, -p - k) in (12TKL The colli- 
sion term (|3-4|) is in the first order of the gradient expansion,^ 1 and the spatial 
point dependence of P a p and Q a p is in the second order. Thus we can adopt the 
coefficients P a p and Q a in the uniform system by neglecting higher order derivative 
terms. 

In the end we can expand the R.H.S. of (|3-3|) by reexpressing by using Dj?, and 

d ^ =9 -r j dFpqk ^ { P ^ (p ' q) [*3r(w*0 - Kl l pT {pqk)] ln|||| 

D 12 (k) ~) 

+ Q^(p,q) [K% L (pqk) - Kf pL (pqk)] ln-^y j , (3-6) 

^ - (ffil (ffii (2^1 W^(P + * + «) ■ (3^) 
K X ^\ P qk) =D 1 a 2 > 21 (X,p)Df> 21 (X, q)D 1 T 2 ' 21 (X, k) , (3-8) 

where we have used the relation D^^—k) = D T 2 L (k). In R.H.S. of (|3-6|) we have 
three types of three gluon interaction processes, that is three transverse (TTT), 
two transverse and one longitudinal (TTL) and one transverse and two longitudinal 
processes (TLL). We have no three longitudinal processes which do not exist for any 
dimensions since Qll = 0. 

For three transverse (TTT) processes, only those terms with Ptt contribute, 
then we apply the symmetry relation (|2-22p . which tells us that Ptt is symmetry 
under any exchange of p, q, k. Thus we can omit the variables in Ptt in the integral, 
Ptt (Pj q) — -Ptt (q,k) = ••• = PTT^pqk). The exchange of variables k o p and 
k f-> r and averaging of the (TTT) part in Eq. (j3 6|) read: 

(TTT) = ^ J dr pqk P TT (pqk) [K T 2 TT (pqk) - K 21 TT (pqk)] In ^r4^ > , 

(3-9) 

because Ptt > and (x — y)ln(x/y) > 0. Here note that entropy production is 
zero for three gluons with collinear momenta due to the form of Ptt- Equality holds 
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when the following relation is satisfied; 

where /3(X) are arbitrary functions for the center of coordinate X. Thus we note 
that the local Bose-Einstein distribution function is realized in equilibrium, 

For (TTL) processes, we have contributions from Plt, Ptl and Qtt terms. By 
exchanging the momentum p f-> k and q -H- k in Plt and Ptl terms, respectively, 
and by using the relation PltO*-, q) = Ptl(Pi k) = Qtt(p, q), we get 

,2 



(TTL) = J d r pqk [K T 2 TL {pqk) - K T l TL (pqk)] 



D T 2 (X,p) n . £>i 2 (X,g) ^ . £>} 2 (X,fc) 

= ^ / dr„ t Q T r(p. q) [JT&itw*) - K? TL ij,qk)] ln > » • (3-12) 

For (TLL) processes, we have contributions from Pll, Qlt and Qtl terms. In 
a similar way to the (TTL) processes, by exchanging the momentum q -H- k and 
p k in Qlt and Qtl terms, respectively, and by using the relation Qx/r(p,k) = 
QtlO*-, q) = Pll (p,q), we get 

(TLL) = ^| dr Mfc P LL (p,q) [KlUpqk) - Kf LT {pqk)] In > . 

(3-13) 

As a result the divergence of the entropy current for the non-Abelian gauge 
theory in the temporal axial gauge is found to be positive definite: 

= (TTT : (1231)) + (TTL : §FT2b ) + (TLL : (ETTBl ) > 0. (3-14) 

Therefore the off-shell processes of gluons are shown to increase the entropy and to 
help the thermalization in the leading order gradient expanded KB equation in TAG. 

There are several comments on the entropy production from off-shell processes 
in order. First, the gluon off-shell processes include -H- 3 and 1 f-> 2 processes. 
These processes are included in KB equation which conserves energy-momentum. It 
is seen by taking quasiparticle approximation (zero spectral width) and remaining 
(keeping?) memory integral as shown in Ref. fl~8|) . These off-shell processes con- 
tribute to the dynamics except for the very late time regions where the distribution 
function change is sufficiently moderate. Secondly, in the leading order of the gradi- 
ent expansion, -H- 3 processes are forbidden by the energy-momentum conservation, 
while they have finite contributions when the memory integral is evaluated explic- 
itly without the gradient expansionPS Thirdly, the 1h2 processes correspond to 
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Landau damping diagrammatically and can be treated as the Vlasov term. In this 
paper we treat them as collision terms without taking the HTL approximation, and 
they are shown to contribute to entropy production. One may wonder why does the 
Vlasov dynamics generate entropy. Entropy production in quantum systems always 
requires coarse graining. The 1 f> 2 processes induce the non-linearity in the equa- 
tion of motion for spectral and statistical functions. As a result, the distribution 
function becomes complex in phase space, and a chaotic behavior will emerge. The 
gradient expansion corresponds to the coarse graining in time, then it induces the 
information loss of the distribution function, namely the entropy production. 

To complete the proof of the H-theorem, we need to solve the gauge dependence 
problem. Around the equilibrium, the controlled gauge dependence was proven in 
Refs. I4"5j) andHHJ). The controlled gauge dependence means that the gauge depen- 
dence always becomes higher order 5T ~ 0(g 2L ) in the coupling expansion under 
the gauge transformation of the truncated effective action (thermodynamic poten- 
tial) r ~ 0(g 2L ~ 2 ), where L is the number of loops. Thus entropy density around 
the thermal equilibrium^ derived by differentiating r (thermodynamic potential) 
with temperature has controlled gauge dependence. Unfortunately there is no proof 
of the controlled gauge dependence under nonequilibrium conditions, and its proof 
is out of the scope of this paper. 

§4. Numerical simulation of the transverse mode 

While there are several problems in numerical simulations as mentioned in the 
Introduction, it is important to investigate the KB dynamics of non-Abelian gauge 
theories numerically. Here we perform the numerical simulation of the transverse 
mode of the gauge theory in the temporal axial gauge in 2+1 dimensions, as a 
first step of non-equilibrium quantum field theoretical investigation of glasma. We 
consider a spatially homogeneous system and solve the dynamical equations numeri- 
cally. In order to trace the time evolution of statistical and spectral functions in the 
Kadanoff-Baym equation, at first we shall use a normal lattice discretization with 
respect to momentum spaceP^ We discretize the spatial length L into 2N S grid 
points with the spatial lattice spacing of a s . Here k 2 is discretized in the following, 

*^±^^t) («) 

where the momenta k n are discretized according to k n = with ni(i = 1,2) £ 
{— N s , ■ ■ ■ ,N S }. By taking k 2 as the above way we can remove most of the lattice 
artifacts We take the following replacement in integrating in the momentum 
space 

where S = 1? is the spatial area. The mesh size N s is chosen to be large enough to 
extract the momentum distribution and to satisfy the conversions of the solutions in 
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each momentum mode. In the numerical analysis we have used a space lattice with 
ma s = 0.3, N s = 30, at/a s = 0.1, Nt = 400 where m is the mass explained in the 
next paragraph, at is a temporal step size and Nt is the number of times for which we 
keep the propagator in memory in order to compute the memory integrals^- Because 
of the exponential damping of unequal-time Green's functions, we can introduce the 
cut-off in the memory integral. 

We adopt some simplifications in actual numerical simulations. We are now 
interested in the dynamics with nonlocal self-energy iI non i. Hence in our simulation 
we replace the tadpole mass shift by constant mass m 2 , which simulates the thermal 
mass, 

#ioc,T -> m 2 - Re7T TiR . (4-3) 

Then we make Green's functions and space-time discretization dimensionless by m. 
The retarded transverse self-energy is given as 

g 2 N f d d p P TT (-p,k + p)(wp + w k+p ) 



4 J (2ir) d WpWp +k [(w p + Wp+k) 2 - (k [) y\ 
w p =v / P 2 Tw 2 (4-5) 

which is necessary to remove logarithmic divergence of the nonlocal self-energy and 
realize the stability of the Green's function for momentum cutoff.®) Furthermore 
we only consider the dynamics of the transverse modes and three transverse (TTT) 
processes. It might be insufficient to take only them in gauge theory, but we will see 
that they play a significant role in the dynamics of KB equation. 

Here we should comment on Ptt function (|2-19ap . This function has ambiguity 
at three points k, p and p + k = where we must evaluate the angle between the 
vector with a zero vector. For example at the limit k — > the function Ptt has 
the possibility to take values in the range < Ptt < 4p 2 . This means that if k 
approaches in the parallel direction with p or p + k, this function takes zero, while 
if in the orthogonal direction, it takes 4p 2 . In these three points, we set Ptt — 
(the least entropy production case; Case I), or we set P TT = 2p 2 ,2(p + k) 2 , and 
2k 2 at k, p and p + k = (Case II), where smoothly changing function is utilized 
in connecting parallel and orthogonal direction around the zero vector. We trace 
the dynamics for both of these cases. In the later discussion, we mainly discuss the 
results in Case II. 

We use the initial condition for Fourier transformed transverse spectral functions 
based on quantum commutation relations 

d x opr(x ,y ;k) 
d x od v op T (x ,y ;k) n =0. (4-6) 







=0 , 


x° 


=y° 








=1 , 


x° 


=y° 








= 0. 


x° 


=y° 





*^ We have checked the stability of our numerical simulations by changing these parameters. 
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We adopt an initial condition of statistical functions which contain the information 
of number distribution functions as 

1 {< + 1/2) 



F T (x ,y ;K 
d x oF T (x ,y°;k) 
d x od yQ F T (x ,y°-X 



=to 







+ 1/2) 



where n£ is initial number distribution functions, and a>£ = \/k 2 + m 2 . 
adopt the so-called "Tsunami" distribution expressed as 



with af/m 2 = 10.4 x , cr 2 /m 2 = 100 x 



^ 6XP 

2 9/2 



1a\ 



ky 

2oi 



(4-7) 
We here 

(4-8) 



' 2tt n 
. mi ^ 



fc t = 9.20 • 2vr/L and TV = 0.25. 



This function has Gaussian peaks around k x = kt and k y = 0, and the widths are 
different for k x and k y directions, as shown in Fig. [31 This might be seen as the toy 
model for gluons in heavy ion collisions. In the following we omit the subscript "T" 
in Green's functions. 




Fig. 3. Initial distribution function n^(rnX° = 0.03). This function is plotted in k x 
plane. 

Now let us discuss the time evolution of the KB equation for the transverse 
modes. We solve f[2-12ah and (12-12bl) with mass m and self-energy, which contains 
three transverse Green's functions. Time evolution of Green's functions is performed 
by referring to and in a similar way as Ref. I44|) . We follow the evolution of n^(X°) 
as a function of oj p defined as 



™ k (^°) + I = x /F(d x od y oF)-(d x oFr , (4-9) 



--X° 
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o> k (A°) = ^/(d x od y0 F)/F 
starting from the initial condition Eq. 



x o =y o =x o 



(4-10) 





Fig. 4. Number distribution function ?ik at mX° — 3.0(upper left), rnX° = 12 (upper right). 
mX° = 30(lower left) and mX° = 90 (lower right). The coupling constant and the number of 
color are g 2 /m = 0.08 and N = 3. 

In Fig. [H we show the time evolution of nk(A°) at 3.0 < mX° < 90 with 
g 2 /m = 0.08 and N = 3. In Fig. [5j we compare the number distribution as a 
function of hi Case I and II. There are several points to be noted from these dis- 
tribution. First, we find that the distribution function converges to the equilibrium 
distribution. At later times, the number distribution function loses its Gaussian 
structure and becomes isotropic in momentum space. Thus we have confirmed that 
off-shell transverse g gg processes can cause kinetic equilibrium. Secondly, the 
equilibrium distribution at low energy is in agreement with the Bose distribution. At 
later times, the temperature and the chemical potential are found to be T ~ 2.8(3.2) 
and n ~ — 1.7(— 1.4) at mX° = 120(174). The chemical potential appears to con- 
verge to zero as expected because the number of gluons is not a conserved variable. 
Simulations starting from different initial conditions with the same energy density 
show the same distributions at late times. Thirdly, high momentum modes are found 
to need longer time to grow up. For the modes with lu > 6, the distribution function 
is still different from the Bose distribution at mX° ~ 120. Finally, we find that 
the results are almost the same for non-zero momentum modes in Cases I and II. 
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The zero momentum mode of n k (X°) in Case I remains the same value since zero 
momentum mode of the self-energy is always zero, but it can evolve in Case II and 
approaches the Bose distribution value at zero momentum. We can conclude that 
the ambiguities at k, p, k + p = does not induce any practical problems. 

Figure [6] shows the time evolution of the entropy density s°/m 2 (|4-lip with 
g 2 /m = 0.08 and 0.06. Here we employ another entropy density expression, 

/d d k 
— p[(l + n k )ln(l + n k )-n k lnn k ], (4-11) 

which can be a criteria of thermalizationpl 1 instead of entropy density (|3-ip to 
reduce numerical cost. The entropy density Eq. (|4-lip is an approximate expres- 
sion of Eq. ()3-ip . and it can be derived from ()3- 1|) in the quasi-particle limit; by 
setting U p -> 0, replacing G 12 = -ipf = 2vr5((A; ) 2 - ojl)(8(-k°) + n k ) and 
G 21 = -ip(l + f) = 2TTd((k ) 2 - U)l)(e(k°) + n k ) as in Ref. and neglecting 
longitudinal part. Entropy density increases monotonically in the whole evolution. 
It increases rapidly at early times when the Gaussian structure collapses, and it 
saturates in the later stage when the kinetic equilibrium is nearly achieved. Mono- 
tonically increasing evolution is consistent with the H-theorem. Thus we confirm 
again that the leading order self-energy with off-shell processes, such as g -H- gg, con- 
tributes to entropy production and equilibration in the numerical simulation without 
gradient expansion in the temporal axial gauge. This entropy production is charac- 
teristic in off-shell KB dynamics, since Boltzmann equation prohibits g <R- gg. This 
property of off-shell effects may help the understanding of the early thermalization 
of dense nonequilibrium gluonic system in the initial stage of the heavy ion collision. 




Fig. 5. Logarithmic plot of number distribution function at late times. Left and right panels 
show the results in Case I and Case II, respectively. In Case I, we show the zero momentum 
mode with the filled circle. 



Before closing, we examine the energy conservation in the present numerical 
simulation. The total energy is given by 

ttot = Ckin + e pot (4-12) 
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e km /L d = \ J^- d (p 2 + m 2 + d x od y o)F(x ,y°;p)\ x0=y0=xo (4-13) 
€ ^ /Ld = \ /(^ /(X °' P) " \ j^en T ^ vac F{X^X^) 

r x ° 

I(X°, p) = / dt' [n p (X°,t'; p)F(t',X°; p) - I7 F (X°, t>; p)p{t',X°; p)] (4-14) 
J o 

as shown in Refs. I23j) . 139]) . l4"0j) . IBTj) . Figure [7] shows the kinetic, potential and total 
energy as a function of time X° for g 2 /m = 0.08 measured from their initial values. 
We show the results relative to the initial total energy density, etot/ m ~ £kin/ m = 
18800 in our lattice spacing. We find that the growth of the kinetic energy is cancelled 
by the decrease of potential energy, and that the energy error is within 0.5 percent. 
This is also true in Case I. 




§5. Summary 



In this paper, first we have reviewed the Kadanoff-Baym (KB) dynamics for 
the non-Abelian gauge theory in the temporal axial gauge with the 2PI effective 
action. As an interaction, we have shown the form of leading order self-energy in 
KB equation. In the hard thermal loop (HTL) approximation, the leading order 
self-energy contributes to a mass shift (Debye mass or the thermal mass). With- 
out the HTL approximation, nonlocal part of the self-energy couples two different 
momentum modes, and have a possibility to produce the width of the spectral func- 
tion and change the number distribution function during the time evolution. By 
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Fig. 7. Evolution of kinetic (solid), potential (dotted) and total energy (bold solid) in units of m 
with g 2 /m = 0.08. 



assuming that the system is spatially uniform and isotropic in the color space, we 
have decomposed the Fourier transformed Green's function into the transverse and 
longitudinal components, and expressed the Schwinger-Dyson (SD) equation in the 
form of the KB equation, i.e. the coupled evolution equation of the spectral and 
statistical functions. 

We have introduced the kinetic entropy current and given a proof of the H- 
theorem for the leading order self-energy. From this proof we find that nonlocal 
part of self-energy can induce entropy production. The entropy is produced via the 
off-shell processes such as g gg, and the transverse part of the entropy produc- 
tion, (TTT) in Eq. (|3-9p . is shown to stop increasing when the distribution function 
converges to the Bose-Einstein distribution. 

Finally we have considered 'practical' time evolution of the transverse modes 
of the KB equation numerically, as a first step of the non-equilibrium quantum 
field theoretical simulation of the non-Abelian gauge theory. We have approximated 
the local part of the self-energy by the constant mass term and the retarded local 
self-energy in vacuum, and the longitudinal part is ignored. We have evaluated 
the time evolution of this practical KB equation and investigated the equilibration 
of distribution functions. During the evolution, the number distribution function 
loses its initial information and approaches the Bose distributions. Off-shell particle 
number changing processes from the leading order self-energy are demonstrated to 
change the number distribution function in numerical simulation without gradient 
expansion. We have confirmed the entropy production from the off-shell processes. 
It is consistent with the proof of H-theorem. 

The Kadanoff-Baym equation has a property that the width of the spectral 
function in the p° space and the memory effects allow off-shell collisions, which vary 
the momentum distribution of particle and induce entropy production. It is natural 
to think that these off-shell effects can contribute more rapid entropy production 
also in the 3+1 dimensions. This would be very important to investigate the entropy 
production in the heavy ion collisions. In order to apply the Kadanoff-Baym equation 
to actual heavy ion collision problems, it is necessary to include higher order effects 
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in the coupling expansion, coupling with the color glass condensate (CGC) which 
appears as the classical field, and 3+1 dimensional numerical calculations. Works in 
these directions are in progress, and will be reported elsewhere. 
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